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Abstract: The effect of chemical reaction on an unsteady magneto hydrodynamic flow past an infinite vertical porous plate 
with variable suction and heat convective mass transfer, where the plate temperature oscillates with the same frequency as 
that of variable suction velocity. The non - linear partial differential equations governing the flow have been solved 
numerically using finite difference method. The flow phenomenon has been characterized with the help of flow parameters 
such as velocity, temperature and concentration profiles for different parameters such as thermal Grashof number (Grashof 
number for heat transfer) (Gt) , solutal Grashof number ( Grashof number for mass transfer) (Gc) , Prandtl number (Pr), 

Schmidt number (5c) , Hartmann number (M), Permeability parameter (K) and Chemical reaction parameter {k r ), 
The velocity, temperature and concentration profiles are shown through graphically and Skin - friction, Nusselt number and 
Schmidt number are shown through tabular forms. 
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Nomenclature 

A Variable suction parameter 

B o Magnetic field component along y — axis 

C Specific heat at constant pressure 

Gr Grashof number 

Gc Modified Grashof number 

g Acceleration of gravity 

K' The permeability of medium 

K The permeability parameter 

M Hartmann number 

Pr Prandtl number 

Sc Schmidt number 

k r Chemical reaction parameter 

D Chemical molecular diffusivity 
T' Temperature of fluid near the plate 

T' Temperature of the fluid far away of the fluid from the plate 

T x Temperature of the fluid at infinity 

C Concentration of the fluid 

C' Concentration of fluid near the plate 

C w Concentration of the fluid far away of the fluid from the plate 

C' m Concentration of the fluid at infinity 

t' Time in x' , y coordinate system 

t Time in dimensionless co - ordinates 

U ' Velocity component in x' — direction 

v' Velocity component in y' — direction 

U Dimensionless velocity component in X — direction 

Nu Nusselt number 
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Sh Sherwood number 
R e Reynold's number 

x', y' Co - ordinate system 

X, y Dimensionless coordinates 

V 0 Mean suction velocity 
Greek symbols 

P Coefficient of volume expansion for heat transfer 

p Coefficient of volume expansion for mass transfer 

£ Smallest positive constant 

K T Thermal conductivity of the fluid 

<7 Electrical conductivity of the fluid 

V Kinematic viscosity 

0 Non - dimensional temperature 

p Density of the fluid 

T Skin - friction 

/U Viscosity, Ns/m 2 

I. Introduction 

Combined heat and mass transfer problems with chemical reaction are of importance in many processes and have, 
therefore, received a considerable amount of attention in recent years. In processes such as drying, evaporation at the surface 
of a water body, energy transfer in a wet cooling tower and the flow in a desert cooler, heat and mass transfer occur 
simultaneously. Possible applications of this type of flow can be found in many industries. For example, in the power 
industry, among the methods of generating electric power is one in which electrical energy is extracted directly from a 
moving conducting fluid. Many practical diffusive operations involve the molecular diffusion of a species in the presence of 
chemical reaction within or at the boundary. There are two types of reactions, homogeneous reaction and heterogeneous 
reaction. A homogeneous reaction is one that occurs uniformly throughout a given phase. The species generation in a 
homogeneous reaction is analogous to internal source of heat generation. In contrast, a heterogeneous reaction takes place in 
a restricted region or within the boundary of a phase. It can therefore be treated as a boundary condition similar to the 
constant heat flux condition in heat transfer. The study of heat and mass transfer with chemical reaction is of great practical 
importance to engineers and scientists because of its almost universal occurrence in many branches of science and 
engineering. The flow of a fluid past a wedge is of fundamental importance since this type of flow constitutes a general and 
wide class of flows in which the free stream velocity is proportional to a power of the length coordinate measured from the 
stagnation point. In many transport processes in nature and in industrial applications in which heat and mass transfer is a 
consequence of buoyancy effects caused by diffusion of heat and chemical species. The study of such processes is useful for 
improving a number of chemical technologies, such as polymer production and food processing. In nature, the presence of 
pure air or water is impossible. Some foreign mass may be present either naturally or mixed with the air or water. The 
present trend in the field of chemical reaction with viscosity analysis is to give a mathematical model for the system to 
predict the reactor performance. A large amount of research work has been reported in this field. In particular, the study of 
chemical reaction, heat and mass transfer is of considerable importance in chemical and hydrometallurgical industries. 
Chemical reaction can be codified as either heterogeneous or homogeneous processes. This depends on whether they occur 
at an interface or as a single phase volume reaction. 

Both free and forced convection boundary layer flows with Soret and Dufour have been addressed by Abreu et al. 
[1]. Afify [2] carried out an analysis to study free convective heat and mass transfer of an incompressible, electrically 
conducting fluid over a stretching sheet in the presence of suction and injection with thermal diffusion and diffusion thermo 
effects. Beg et al. [3] have also analyzed the chemical reaction rate effects on steady buoyancy - driven dissipative 
micropolar free convective heat and mass transfer in a Darcian porous regime. Beg et al. [4] have also studied numerically 
both viscous heating and Joule (Ohmic) dissipation effects on transient Hartmann - Couette convective flow in a Darcian 
porous medium channel also including Hall current and ion - slip effects. Chaudhary et al. [5] studied the effect of free 
convection effects on magnetohydrodynamic flow past an infinite vertical accelerated plate embedded in porous media with 
constant heat flux by using Laplace transform technique for finding the analytical solutions. Chin et al. [6] obtained 
numerical results for the steady mixed convection boundary layer flow over a vertical impermeable surface embedded in a 
porous medium when the viscosity of the fluid varies inversely as a linear function of the temperature. Das and Mitra [7] 
discussed the unsteady mixed convective MHD flow and mass transfer past an accelerated infinite vertical plate with 
suction. Recently, Das and his co - workers [8] analyzed the effect of mass transfer on MHD flow and heat transfer 
past a vertical porous plate through a porous medium under oscillatory suction and heat source. Das et al. [9] 
investigated numerically the unsteady free convective MHD flow past an accelerated vertical plate with suction and heat 
flux. Das and his associates [10] estimated the mass transfer effects on unsteady flow past an accelerated vertical 
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porous plate with suction employing finite difference analysis. Gireesh kumar et al. [11] investigated effects of 
chemical reaction and mass transfer on MHD unsteady free convection flow past an infinite vertical plate with 
constant suction and heat sink. 

Hayat et al. [12] analyzed a mathematical model in order to study the heat and mass transfer characteristics in 
mixed convection boundary layer flow about a linearly stretching vertical surface in a porous medium filled with a visco - 
elastic fluid, by taking into account the diffusion thermo and thermal diffusion effects. Ibrahim et al. [13] studied the effects 
of chemical reaction and radiation absorption on transient hydromagnetic natural convection flow with wall transpiration and 
heat source. The unsteady free convective MHD flow with heat transfer past a semi - infinite vertical porous moving 
plate with variable suction has been studied by Kim [14]. Effects of chemical reaction on free convection flow of a polar 
fluid through a porous medium in the presence of internal heat generation are examined by Patil and Kulkarni [15]. Several 
studies have also described thermal radiation effects on convection flows in porous media. MHD convective flow of a 
micropolar fluid past a continuously moving vertical porous plate in the presence of heat generation/absorption was studied 
by Rahman and Sattar [16]. Sarangi and Jose [17] studied the unsteady free convective MHD flow and mass transfer 
past a vertical porous plate with variable temperature. The effect of temperature dependent viscosity and thermal 
conductivity on unsteady MHD convective heat transfer past a semi - infinite vertical porous plate has studied Seddek and 
Salama [18]. Seddeek et al. [19] more recently reported on the effect of chemical reaction and variable viscosity on 
hydromagnetic mixed convection heat and mass transfer for Hiemenz flow through a Darcian porous media in the presence 
of radiation and magnetic field. Seethamahalakshmi et al. [20] discussed the effects of Soret on an unsteady 
magnetohydrodynamic free convection flow past an infinite vertical porous plate with variable suction. The unsteady free 
convection flow of a viscous incompressible fluid past an infinite vertical plate with constant heat flux is considered on 
taking into account viscous dissipative heat, under the influence of a transverse magnetic field studied by Srihari et al [21]. 

The object of the present paper is to analyze the chemical reaction effect on an unsteady magneto hydrodynamic 
free convection flow past an infinite vertical plate by taking variable suction into account. The governing equations are 
transformed by using unsteady similarity transformation and the resultant dimensionless equations are solved by using the 
finite difference method. The effects of various governing parameters on the velocity, temperature, concentration, skin - 
friction coefficient, Nusselt number and Sherwood number are shown in figures and tables and discussed in detail. From 
computational point of view it is identified and proved beyond all doubts that the finite difference method is more 
economical in arriving at the solution and the results obtained are good agreement with the results of Seethamahalakshmi et 
al. [20] in some special cases. 

II. Mathematical Analysis 

Consider the unsteady magneto hydrodynamic flow past an infinite vertical porous plate with variable suction and 
heat convective mass transfer, where the plate temperature oscillates with the same frequency as that of variable suction 
velocity in presence of chemical reaction. 
We made the following assumptions. 

1 . In Cartesian coordinate system, let X — axis is taken to be along the plate and the y ' — axis normal to the plate. Since 
the plate is considered infinite in X — direction, hence all physical quantities will be independent of X — direction. 

2. The wall is maintained at constant temperature (7"') and concentration (C') higher than the ambient temperature 
(7^) and concentration (C^) respectively. 

3. Let the components of velocity along x' and y' axes be u' and v' which are chosen in the upward direction along the 
plate and normal to the plate respectively. 

4. A uniform magnetic field of magnitude B g is applied normal to the plate. The transverse applied magnetic field and 

magnetic Reynold's number are assumed to be very small, so that the induced magnetic field is negligible. 

5. The polarization effects are assumed to be negligible and hence the electric field is also negligible. 

6. The homogeneous chemical reaction of first order with rate constant between the diffusing species and the fluid is 
assumed. 

7. The concentration of the diffusing species in the binary mixture is assumed to be very small in comparison with the 
other chemical species, which are present and hence Soret and Dufour effects are negligible. 

8. The Hall Effect of magneto hydrodynamics and magnetic dissipation (Joule heating of the fluid) are neglected. 



Hence the governing equations of the problem are: 
dp' { d(p'u') | d(p'u') 
8t' 
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Here, the status of an equation of state is that of equation p = constant. This means that the density variations 

produced by the pressure, temperature and concentration variations are sufficiently small to be unimportant. Variations of all 
fluid properties other than the variations of density except in so far as they give rise to a body force, are ignored completely 
(Boussinesq approximation). All the physical variables are functions of y 'and f'only as the plate are infinite. It is also 

assumed that the variation of expansion coefficient is negligibly small and the pressure and influence of the pressure on the 
density are negligible. Within the framework of above assumptions the governing equations reduce to Continuity Equation: 



8v' 



= 0 



2 I 



dy' 

Momentum Equation: 

dt dy dy 

Energy Equation: 

dT' ,dV _ K T d 2 T' 
ltf +V ly~~fJC p ~dy 

Species Diffusion Equation: 



V 



P K' 



i2 



dt' dy' 



dy 
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(5) 



(6) 



(7) 



(8) 



And the corresponding boundary conditions are 
t'<0:u' = 0, T' = Z, C' = C' m for all y' 



t'>0:- 



u' — 0, T' — T' w — 1 + se imY , C' = C' W aty' = 0 
w'->0, T'^T^, C'^C' X as /->oo 



(9) 



From the continuity equation, it can be seen that v' is either a constant or a function of time. So, assuming suction 
velocity to be oscillatory about a non - zero constant mean, one can write V = — V a (l + sAe 1 ' 0 ''') 

(10) 

Where £, A are small such that £"A<<1. The negative sign indicates that the suction velocity is directed towards the 
plate. In order to write the governing equations and the boundary conditions in dimensional following non - dimensional 
quantities are introduced. 
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Hence, using the above non - dimensional quantities, the equations (6) - (9) in the non - dimensional form can be written as 



1 du 



(l + eAe ^ \Su = du + (Gr)0 + (Gc)C - 
4 dt V 'dy dy 2 
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i ee 



\d0 1 d 2 0 



4 dt dy Pr dy 2 



l - d -^-{l + sAe-')°— = -^-k r C 
4 dt 'dy Sc dy 2 

And the corresponding boundary conditions are 
t < 0 : u = 0, 0 = 0, C = 0 for all y 

r u=0, 0 = l + se" a ,C = l at y = 0 



\dC 1 d 2 C 



(13) 
(14) 



t>0:< 



•0, 6>^0, C^O as y 



• 00 



(15) 



All the physical parameters are defined in the nomenclature. 

It is now important to calculate the physical quantities of primary interest, which are the local wall shear stress, the 
local surface heat and mass flux. Given the velocity field in the boundary layer, we can now calculate the local wall shear 
stress (i.e., skin - friction) is given by and in dimensionless form, we obtain Knowing the temperature field, it is interesting 
to study the effect of the free convection and radiation on the rate of heat transfer. This is given by which is written in 
dimensionless form as 



du 

dy 



y'=0 



du 
~dy~ 



(16) 



v=0 



The dimensionless local surface heat flux (i.e., Nusselt number) is obtained as 



N u (x') = 



x 



dT 



(K-K) dy' 



Then Nu = 



de_ 

dy 



(17) 



v=0 



J/=o 

The definition of the local mass flux and the local Sherwood number are respectively given by with the help of these 
equations, one can write 



S h (x') 



X 



dc 



(C' w -C'Jdy' 



Then Sh = 



y'=0 



S h (x') 
R 



dc_ 

dy 



(18) 



v=() 



n V o X 

Where K e = the Reynold is's number. 

V 

III. Method of solution 

Equations (12) - (14) are coupled non - linear partial differential equations and are solved by using initial and 
boundary conditions (15). However, exact or approximate solutions are not possible for this set of equations. And hence we 
solve these equations by an implicit finite difference method of Crank - Nicolson type for a numerical solution. The 
equivalent finite difference scheme of equations (12) - (14) is as follows: 
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Here the suffix i corresponds to y and j corresponds to t . Also 5 = 1 + sAe'" 1 , At = t . +J — t . and Ay = y M — y i 
The complete solution of the discrete equations (19) - (21) proceeds as follows: 

Step - (1): Knowing the values of C , 9 and u at a time t = j , calculate C and 9 at a time t = j + 1 using equations 
(20) and (21) and solving tri - diagonal linear system of equations. 

Step - (2): Knowing the values of 9 and C at time? = j , solve the equation (19) (via tri - diagonal matrix inversion) to 
obtain u at a time t = j + 1 . 

We can repeat steps (1) and (2) to proceed from t — 0 to the desired time value. 

The implicit Crank - Nicolson method is a second order method [pi/St 2 )) in time and has no restrictions on space 
and time - steps, Ay and A?, i.e. ,the method is unconditionally stable. The finite differences scheme used, involves the 
values of the function at the six grid points. A linear combination of the "future" points is equal to another linear 
combination of the "present" points. To find the future values of the function, one must solve a system of linear equations, 
whose matrix has a tri - diagonal form. 

The computations were carried out for Gr — 1.0, Gc — 1.0, Pr =0.71 (Air), Sc - 0.22 (Hydrogen), M = 1.0, 

7t 

# = 1.0, A = 1.0, £ = 0.001, Q)t = — ,k r = 1.0, t =1.0 and Ay =0.1, At = 0.001 and the procedure is repeated till y = 

4. In order to check the accuracy of numerical results, the present study is compared with the available theoretical solution of 
Seethamahalakshmi et al. [20] and they are found to be in good agreement. 

IV. Results and Discussions 

In the preceding sections, the problem of chemical reaction effects on an unsteady MHD free convection flow past 
an infinite vertical porous plate with constant suction was formulated and the dimensionless governing equations were 
solved by means of a finite difference method. In the present study we adopted the following default parameter values of 

finite difference computations: Gr= 1.0, Gc= 1.0, Pr=0.71 (Air), Sc = 0.22 (Hydrogen), M = 1.0, K = 1.0, A = 1.0, 
S = 0.001, COt = — , k r =1.0 and t = 1.0. All graphs therefore correspond to these values unless specifically indicated on 
the appropriate graph. 

The influence of thermal Grashof number, Gr , on the velocity is shown in Figure (1). The thermal Grashof 
number signifies the relative effect of the thermal buoyancy force to the viscous hydrodynamic force. The flow is accelerated 
due to the enhancement in buoyancy force corresponding to an increase in the thermal Grashof number. The positive values 
of Gr correspond to cooling of the plate by natural convection. Heat is therefore conducted away from the vertical plate into 
the fluid which increases the temperature and thereby enhances the buoyancy force. In addition, it is seen that the peak 
values of the velocity increases rapidly near the plate as thermal Grashof number increases and then decays smoothly to the 
free stream velocity. Figure (2) presents typical velocity profiles in the boundary layer for various values of the solutal 
Grashof number, Gc . The solutal Grashof number Gc defines the ratio of the species buoyancy force to the viscous 
hydrodynamic force. It is noticed that the velocity increases with increasing values of the solutal Grashof number. 

Figures (3) and (8) show the behavior of velocity and temperature profiles for different values Prandtl number. The 
numerical results show that the effect of increasing values of Prandtl number results in a decreasing velocity. It is observed 
that an increase in the Prandtl number results a decrease of the thermal boundary layer thickness and in general lower 

average temperature within the boundary layer. The reason is that smaller values of Pr are equivalent to increase in the 
thermal conductivity of the fluid and therefore, heat is able to diffuse away from the heated surface more rapidly for higher 

values of Pr . Hence in the case of smaller Prandtl number as the thermal boundary later is thicker and the rate of heat 
transfer is reduced. The effect of Schmidt number Sc on the velocity and concentration are shown in figures (4) and (9). As 
the Schmidt number increases, the velocity and concentration decreases. This causes the concentration buoyancy effects to 
decrease yielding a reduction in the fluid velocity. Reductions in the velocity and concentration distributions are 
accompanied by simultaneous reductions in the velocity and concentration boundary layers. 

The effect of Hartmann number, M on the velocity is shown in figure (5). The velocity decreases with an increase 
in the Hartmann number. It is because that the application of transverse magnetic field will result a resistive type force 
(Lorentz force) similar to drag force which tends to resist the fluid flow and thus reducing its velocity. Also, the boundary 
layer thickness decreases with an increase in the Hartmann number. Figure (6) shows the effect of the permeability of the 
porous medium parameter, K on the velocity distribution. As shown, the velocity is increasing with the increasing 
dimensionless porous medium parameter. The effect of the dimensionless porous medium K becomes smaller as K 
increase. Physically, this result can be achieved when the holes of the porous medium may be neglected. Figures (7) and 

(10), illustrate the behavior velocity and concentration for different values of chemical reaction parameter, k r . It is observed 
that an increase in leads to a decrease in both the values of velocity and concentration. A distinct velocity escalation occurs 
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near the wall after which profiles decay smoothly to the stationary value in free stream. Chemical reaction therefore boosts 
momentum transfer, i.e., accelerates the flow. 

The profiles for skin - friction (t) due to velocity under the effects of Grashof number (Gr) , Modified Grashof 

number (Gc) , Prandtl number (Pr), Schmidt number (Sc) , Hartmann number (M), Permeability parameter (K) and 

Chemical reaction parameter (k r ) are presented in the table - 1 respectively. We observe from this table - 1, the skin - 

friction (r) increases under the effects of Grashof number (Gr), Modified Grashof number (Gc) and Permeability 

parameter (K) . And decreases under the effects of Prandtl number(Pr), Schmidt number(Sc) , Hartmann number (M) 

and Chemical reaction parameter (k r ) 



Table - 1: Skin - friction results [T ) for the values of Gr, Gc, Pr, Sc, M, K and k r 



Gr 


Gc 


Pr 


Sc 


M 


K 


K 


T 


1.0 


1.0 


0.71 


0.22 


1.0 


1.0 


1.0 


5.8411 


2.0 


1.0 


0.71 


0.22 


1.0 


1.0 


1.0 


8.9829 


1.0 


2.0 


0.71 


0.22 


1.0 


1.0 


1.0 


8.5404 


1.0 


1.0 


7.00 


0.22 


1.0 


1.0 


1.0 


4.4432 


1.0 


1.0 


0.71 


0.30 


1.0 


1.0 


1.0 


5.7113 


1.0 


1.0 


0.71 


0.22 


2.0 


1.0 


1.0 


4.5870 


1.0 


1.0 


0.71 


0.22 


1.0 


2.0 


1.0 


6.8033 


1.0 


1.0 


0.71 


0.22 


1.0 


1.0 


2.0 


5.5172 



Table - 2: Rate of heat transfer ( Nu ) values for different values of Pr and Rate of mass transfer ( Sh ) values for different 



values of Sc and k r 



Pr 


Nu 


Sc 


K 


Sh 


0.71 


7.7337 


0.22 


1.0 


6.7401 


7.00 


4.9322 


0.30 


1.0 


6.4490 


0.22 


2.0 


6.0074 



The profiles for Nusselt number (Nu) due to temperature profile under the effect of Prandtl number (Pr)is 
presented in the table - 2. From this table we observe that, the Nusselt number (Nu) due to temperature profiles decreases 
under the effect of Prandtl number (Pr). The profiles for Sherwood number (Sh) due to concentration profiles under the 
effect of Schmidt number (Sc) and Chemical reaction parameter (k r ) are presented in the table - 2. We see from this table 
the Sherwood number (Sh) due to concentration profiles decreases under the effects of Schmidt number (Sc) and 
Chemical reaction parameter (k r ) . 

Table - 3: Comparison of present Skin - friction results ( r ) with the Skin - friction results ( T ) obtained by 
Seethamahalakshmi et al. [20] for different values of Gr, Gc, Pr, Sc, M and K 



Gr 


Gc 


Pr 


Sc 


M 


K 


T 


T* 


1.0 


1.0 


0.71 


0.22 


1.0 


1.0 


4.5517 


4.5508 


2.0 


1.0 


0.71 


0.22 


1.0 


1.0 


6.2941 


6.2931 


1.0 


2.0 


0.71 


0.22 


1.0 


1.0 


6.5403 


6.5397 


1.0 


1.0 


7.00 


0.22 


1.0 


1.0 


3.9954 


3.9922 


1.0 


1.0 


0.71 


0.30 


1.0 


1.0 


4.2225 


4.2216 


1.0 


1.0 


0.71 


0.22 


2.0 


1.0 


3.8467 


3.8449 


1.0 


1.0 


0.71 


0.22 


1.0 


2.0 


4.9277 


4.9260 



In order to ascertain the accuracy of the numerical results, the present results are compared with the previous results 
of Seethamahalakshmi et al. [20] for Gr= Gc= 1.0, Pr = 0.71, Sc = 0.22, M = 1.0 and K= 1.0 in table - 3. They are 
found to be in an excellent agreement. 
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V. Conclusions 

In this paper, the chemical reaction effects on an unsteady MHD free convection flow past an infinite vertical 
porous plate with variable suction was considered. The non - dimensional governing equations are solved with the help of 
finite difference method. The conclusions of the study are as follows: 

1. The velocity increases with the increase in thermal Grashof number, solutal Grashof number and Permeability 
parameter. 

2. The velocity decreases with an increase in the Hartmann number. 

3. An increase in the Prandtl number decreases the velocity and temperature. 

4. The velocity as well as concentration decreases with an increase in the Schmidt number. 

5. The velocity as well as concentration decreases with an increase in the chemical reaction parameter. 

6. The skin - friction increases with the increase in thermal Grashof number, solutal Grashof number and Permeability 
parameter. 

7. The skin - friction decreases with an increase in the Hartmann number, Prandtl number, Schmidt number and Chemical 
reaction parameter. 

8. An increase in the Prandtl number decreases the Nusselt number. 

9. The Sherwood number decreases with an increase in the Schmidt number and Chemical reaction parameter. 

10. On comparing the skin - friction (r) results with the skin - friction (T ) results of Seethamahalakshmi et al. [20] it can 
be seen that they agree very well. 
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Figure 1 . Velocity profiles for different values of Gr 
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Figure 4. Velocity profiles for different values of Sc 
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Figure 5. Velocity profiles for different values of M 



1 




o \ y 2 3 



Figure 6. Velocity profiles for different values of K 
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Figure 7. Velocity profiles for different values of k r 
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Figure 8. Temperature profiles for different values of Pr 
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Figure 9. Concentration profiles for different values of Sc 
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Figure 10. Concentration profiles for different values of k r 
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